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Abstract 

We study the role of fluctuations in single molecule experimental measure- 
ments of force-extension (/ — Q) curves. We use the Worm Like Chain (WLC) 
model to bring out the connection between the Helmholtz ensemble charac- 
terized by the Free Energy {F{C,)) and the Gibbs ensemble characterized by 
the Free Energy {G{f )). We consider the rigid rod limit of the WLC model 
as an instructive special case to bring out the issue of ensemble inequivalence. 
We point out the need for taking into account the free energy of transition 
when one goes from one ensemble to another. We also comment on the "phase 
transition" noticed in an isometric setup for semiflexible polymers and pro- 
pose a realization of its thermodynamic limit. We present general arguments 
which rule out non-monotonic force-extension curves in some ensembles and 

note that these do not apply to the isometric ensemble. 
PACS numbers: 87.15.-v,05.40.-a,36.20.-r 
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Introduction: 

In the past, experiments on polymers were confined to studying their bulk properties, 
which involved probing large numbers of molecules [1]. The results of these experiments 
could be analyzed by using the traditional tools of thermodynamics. In recent years, how- 
ever, researchers have been successful in micromanipulating single biological molecules such 
as DNA, proteins and RNA to probe their elastic properties [2]. Such studies serve a twofold 
role. On one hand, they shed light on mechanical properties of semiflexible polymers, which 
are of clear biological importance in processes such as gene regulation and transcription 
[3-5]. On the other hand, they provide physicists with a concrete testing ground for under- 
standing some of the fundamental ideas of statistical mechanics. In statistical mechanics, 
an isometric setup would be described by the Helmholtz free energy, whereas an isotensional 
setup would be described by the Gibbs free energy [6]. In the thermodynamic limit these 
two descriptions agree, but semiflexible polymers (those with contour lengths comparable to 
their persistence lengths) are not at the thermodynamic limit. Experimentally, both isomet- 
ric and isotensional ensembles are realizable. Typically the polymer molecule is suspended 
(in a suitable medium) between a translation stage and a force sensor. The force sensor 
could be realized by using an Atomic Force Microscope (AFM) cantilever or by optical or 
magnetic forces. As noted by Kreuzer and Payne [7], an isometric setup can be realized 
using a stiff trap and an isotensional setup by using a soft trap. In a more sophisticated 
version, an electronic feedback circuit is used to control the force (or the extension) and 
one measures the fluctuations in the extension (or the force) [8]. Here we will focus on 
the role of fluctuations in single molecule experiments. In order to correctly interpret such 
experiments one needs to understand the effect of fluctuations on the measured quantities. 
For instance, it turns out, that an experiment in which the ends of a polymer molecule are 
fixed (isometric) and the force fluctuates yields a different result from one in which the force 
between the ends is held fixed (isotensional) and the end-to-end distance fluctuates [8,9]. 
This difference can be traced to large fluctuations about the mean value of the force or 
the extension, depending on the experimental setup. These fluctuations vanish only in the 
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thermodynamic limit of very long polymers. 

Here we use the Worm Like Chain (WLC) model [10-12] to study the inequivalence of 
ensembles due to finite size effects. The WLC model has been very successful in achieving 
quantitative agreement with experimentally measured force-extension curves [2,13]. The 
paper is organized as follows. In Sec. // we discuss the Helmholtz and Gibbs ensembles. 
In Sec. /// we consider the rigid rod limit which forcefully brings out the main issues dealt 
with in this paper. In Sec. IV we draw attention to the importance of taking into account 
the free energy of transition in going from one ensemble to another. In Sec. V we discuss the 
thermodynamic limit of a "phase transition" recently seen in semiflexible polymers. Finally, 
we end the paper with a discussion in Sec. VI. 

II Helmholtz and Gibbs Ensembles: 

Consider an idealized experiment in which one end of a molecule is held fixed at {xq, i/q, 0) 
and the other end is attached to a dielectric bead confined to a harmonic optical trap 
described by the potential 

V(x, y, z) = A^i-^^^l±^^y^ + C^l^. (1) 

with (xo, Uo, Zq) defining the center of the trap. Consider A to be small so that the bead is free 
to move in the plane z = Zq. For a polymer of contour length L and persistence length Lp 
it is convenient to introduce the following dimensionless variables: C — 
and / = 1^ where F is an applied stretching force and ksT is the thermal energy at a 
temperature T. Consider P[Qd(, the number of configurations (counted with Boltzmann 
weight) [1 1] for a polymer of length L starting from the origin and ending anywhere on the 
X — y plane in an interval d( of (. The free energy defined by J^{C) — —^lnP{() is the 
Helmholtz free energy. The partition function [14] for the combined system consisting of the 
polymer molecule and the trap is given by: 

^(Co, /3) = # r dCe-^^^'^e-^"^ (2) 

V ZTT J-oo 

where — 01?. By tuning the longitudinal stiffness C we can realize the two limiting cases. 
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Helmholtz : In the limit of a stiff trap {C — > oo), the Gaussian factor pertaining to the trap 
approaches a delta function and one gets 

Z(C,/5) = e-^^K) (3) 

Here we have switched notation to write ( in place of Co- Thus a stiff trap realizes the 
Helmholtz ensemble by constraining fluctuations in the ( coordinate. To extend the molecule 
from C to C + c^C 0116 needs to apply a force < / >= ^ in order to compensate for the 
change of entropy. Plotting < / > versus C we find the (< / >, C) force-extension relation. 

Gibbs : In the opposite limit of a soft trap (C* — > and (q ^ oo such that C(o — f3f 
remains finite), one gets [14] 

Z(/,/3)= pdCe-^^^^^e/^/^ (4) 

Thus a soft trap permits fiuctuations in the C, coordinate but constrains the force fiuctua- 
tions [7] and thus realizes the Gibbs ensemble. Z{f) is the generating function for the C 
distribution. Defining the Gibbs free energy G(/) = —pinZ{f) we can work out the mean 
extension < ( >— — ^ and the (< ( >, /) force-extension relation. 

Notice that Z{f) is the Laplace transform of Z{(). In the thermodynamic limit of long 
polymers (/? —>■ oo) the Laplace transform integral Eq.(4) is dominated by the saddle point 
value and therefore ^(C) and G{f) are related by a Legendre transform: 

HO-G{f) + K. (5) 

For finite (3 i.e. for a polymer of finite extent, the saddle point approximation no longer 
holds true and fiuctuations about the saddle point value of the free energy become important. 
Thus one finds that J-'{C) and G{f ) are not Legendre transforms of each other. We notice 
that this difference between the Legendre transform [Eq. (5)] and the Laplace transform 
[Eq. (4)] is the mathematical origin of the finite size fiuctuation effects described in this 
paper. These fiuctuations are of thermal origin and can ultimately be traced to collisions 
of the polymer molecule with the molecules of the suspending medium. In this section we 
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have recovered the results of [7] . We have also gone beyond [7] and traced the origin of the 
difference between the two ensembles to the difference between the Legendre and Laplace 
transforms. 

Ill Rigid Rod Limit: An Instructive Extreme Case 

We noticed in the last section that because of fluctuation effects the Helmholtz and Gibbs 
free energies are not related by a Legendre Transform. An important consequence of this is 
that the (F, < z >) relation is different from the (< F >,z) relation. In other words, the 
force-extension curves plotted in the two ensembles are distinct due to flnite size fluctuation 
effects. Fluctuations dominate at finite /5 and disappear in the thermodynamic limit {(3 —>■ 
oo) of flexible polymers. We bring out the ensemble dependence of the force-extension 
relations explicitly and most dramatically in the limiting case of a very stiff polymer {/3 ^ 0). 

In this extremely rigid limit [15], the end of the polymer is uniformly distributed over the 
sphere of directions. In the Helmholtz ensemble we thus have P{z)dz — || for —L < z < L 
and P{z) — otherwise. Since the free energy is constant in the range —L < z < L and 
diverges otherwise, we find that the average force < F > vanishes for \z\ < L and diverges 
for \z\ — L (See Fig. la). In the Gibbs ensemble we find by standard manipulations that 

<z>^ [LCothFL - h (6) 

F 

which differs from the < F > —z relation determined above in the Helmholtz ensemble (See 
Fig. 16). The theoretical analysis of the ensemble dependence of the force-extension relation 
based on the rigid rod limit is a new result of this paper. 

Thus an experimenter making force-extension measurements on, for instance, Actin fila- 
ments [16] , would find that a measurement in which the force is controlled and the end-to-end 
distance is measured leads to a different force-extension curve from a measurement in which 
the end-to-end separation is controlled and the force is measured. A theorist interpreting 
the curves also needs to keep in mind whether the curves are obtained in a constant-force 
setup or a constant-extension setup since a proper interpretation of the curves requires a 
knowledge of the ensemble used in the measurement. 

5 




0.2 0.4 0.6 0.8 7 1 0^2 0^4 ol 0^8 

^ <Z> 

FIGURES 



FIG. 1. Force-Extension Curve in the Helmholtz (Fig. la) and Gibbs (Fig. lb) ensembles for 
(3 = 0. We have set L = 1. 

IV Inequivalence Of Ensembles and the Free Energy Of Transition 

In an isotensional setup one controls the force and one measures the mean extension and 
plots it as a function of force. In an isometric setup the roles of extension and force are 
interchanged. In both setups zero force corresponds to zero extension (^ = 0) and a large 
force corresponds to maximal extension (i.e. ( = 1). Imagine going from zero force to a 
large force via the isotensional setup and returning from maximal extension to zero extension 
via the isometric setup. Since the "equation of state" depends on the chosen ensemble, in 
general there will be two distinct curves describing the extension in one ensemble followed 
by contraction in the other [17]. In such a situation there could be a net area enclosed in 
the Force-Extension plane. This poses a puzzle because it appears that a cyclic process can 
extract work from the system. This puzzle is easily resolved. In completing the cycle and 
returning to the initial state one is in fact changing ensembles twice at the two end points. 
These correspond to finite free energy changes which need to be taken into account. 

Let (Ci, /i) and (^2, /2) be two points which lie on both isometric and isotensional force- 
extension curves. In our example (Ci,/i) = (0,0) and {(2,12) = (l,cxo). Let us suppose 
that we go from (Ci, /i) to {(,2-1 f 2) in the isotensional ensemble and return via the isometric 
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ensemble. We find that 

Gu.)-G(U)=j';^'^d!=-j';^<c>d!. (7) 

Similarly, 

rC,2 Qp K2 

F(C.) - = l^^cK = l^<f> dC (8) 

The area enclosed between the two curves is given by 

W= fd<C>- <f>dc 

which can be rewritten as 

W = d{f <(>)- <C>df- f^' <f>dC (9) 

^Ci Jfi ki 

= /2C2 - /iCi + ^(72) - - F(C2) + F(Ci) 

-[K + G{f)-F{C)]l = [F{C)-F{0]l 

= AF*^(C)|? 

where we define AF*^(C) as the free energy of transition. AF*^{Q is the difference 
between F{() = fC + G{f), the Legendre transform of G{f) and the Helmholtz free energy 
F{Q. Since these are not equal (except in the limit of long polymers) the free energy of 
transition between ensembles must be considered in order that the total free energy change 
in a cyclic process vanishes. 

In order to understand this issue more explicitly we consider corrections to the saddle 
point approximation which is valid in the long polymer limit. Let us expand (f){Q — J-{C) — 
fC, the argument of the exponential appearing on the right hand side of Eq. (4) around 
the saddle point value C = C* (which dominates the integral in the long polymer limit) and 
retain terms upto second order in the fluctuations about the saddle point value: 

If we plug in this expansion in Eq. (4) and identify Z(f,(3) with e"^'^*-'^^ we arrive at the 
following equation: 
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G{f) = [HQ-fQ + 7^ln:F"{Q 



(10) 



The free energy due to fluctuations around the saddle point value is i^lnj^ (C*)- Notice 
that in the long polymer limit of /3 ^ oo, this term vanishes. For finite j3, this nonzero 
contribution to the free energy accounts for the transition between the constant extension 
ensemble and the constant force ensemble. In going from a soft trap to a stiff trap work 
is done on the bead by the trap. Similarly in going from a stiff trap to a soft trap work is 
extracted from the bead by the trap. The net work done is the difference between the work 
done at the two ends of the force-extension curves in switching ensembles. This net work 
exactly cancels out the nonzero area enclosed in the force-extension plane. 

V "First Order Phcise Transition" and the Thermodynamic Limit in Semi- 
flexible Polymers 

In Sec. // we considered the bead to be in a potential which was soft in the x and y 
directions. Let us now consider what happens when the trap is stiff in all three directions {A 
as well as C in Eq. (1) are large) and the vector position of the bead is constrained to be at 
(xo, 2/0, zq). Let Q{r) be the number of polymer configurations which start at the origin and 
end in the volume element centered at r [11,18]. Q{f) is related to P{z) via the equation 



which in words, means that P{z) is obtained by integrating Q{r) over a plane of constant z 
[11]. The distribution Q{r) was studied in [18] where it was noticed that in an intermediate 
range (around 3.8) of (3 the free energy A{r) = ^logQif) had multiple minima. For 
a fixed contour length as one varies (3 by tuning the persistence length Lp one finds a 
competition between fiexible and rigid phases of the polymer for intermediate values of (5. 
Thus the polymer undergoes a fiexible to rigid "first order phase transition" via a two peaked 
profile ol Q{r). This leads to a curious force-extension relation. As one pulls the bead, the 
restoring force at first increases, then decreases to zero and then goes negative and becomes 
a destabilizing force. The molecule is unstable and goes to a new extended state. 




(11) 
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The words "first order phase transition" above were in quotes as a finite system does 
not show phase transitions. If one takes the thermodynamic hmit by taking the length of 
the polymer to infinity {(3 — > oo) one loses the multiple minima structure which is present 
only in a small range of /3 around 3.8. Is it possible for this phase transition to survive the 
thermodynamic limit? As we will see below, this is indeed possible provided one takes the 
replica thermodynamic limit. We take N replicas of the molecule with fixed (3 and let N 
tend to oo. Consider N identical polymers with P — fixed, their two ends anchored to 
fiat surfaces 5"! and S2 (See Fig. 2). One could realize the above arrangement by using (i) 
two planar arrays of optical traps or (ii) by introducing suitably synthesized supramolecular 
lamellar structures. The anchoring is such that the tangent vectors to the molecule at 
the fixed ends are free to swivel. If one apphes a force F to pull 5"! and S2 apart the N 
molecules are also stretched. We consider the molecules to be well separated so that they can 
be regarded as independent. One could now look at the mean force F needed to maintain 
the separation r. It is easily seen that the force is proportional to A'^ and also the mean 
square fluctuation < (AF)^ > in the force is proportional to A^. This is because the mean 
force and its variance are respectively the first and second derivatives of the free energy, 
which being an extensive quantity is proportional to A". It follows that as A^ goes to infinity, 
the fluctuations (AF/F) in F die out as 1/\/N. We can now regard the mean force F as 
a control parameter (i.e. consider a constant F ensemble) and observe that if we tune the 
applied force F, there is a discontinuous change in the separation r between the two sets of 
optical traps signalling a flrst order phase transition with the inter-trap separation r as the 
order parameter. Thus, in the replica way of taking the thermodynamic limit the double 
humped form of the distribution function Qlr) results in a true first order phase transition. 
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FIG. 2. Schematic Experimental Design for Replica Thermodynamic Limit, shown above for 
iV = 4. 

VI Discussion: 

In this paper we point out the importance of considering the free energy of transition 
in going between the Helmholtz and the Gibbs ensembles in the context of single molecule 
force-extension measurements. We also study two distinct ways of taking the thermodynamic 
limit (i) by letting the length of the polymer tend to infinity (i. e. /3 — > oo) and {ii) by 
considering replicas. In particular we notice that the flexible to rigid transition mentioned 
in Sec. V for a single semifiexible polymer survives the replica thermodynamic limit. In 
contrast, this feature disappears in the usual thermodynamic limit of /3 — > oo. 

The non-monotonic behavior of Qir) is an intriguing feature of molecular elasticity. It 
was noticed in recent simulations [18] and subsequently in a semi-analytical treatment [11]. 
It was also commented on in a recent paper [19]. This double humped form of Q{r) has 
gained renewed interest in the context of cyclization of polymers and its significance to gene 
regulation [20,21]. 

We note that this remarkable qualitatively distinctive feature is special to the isometric 
ensemble (fixed vector end-to-end separation) and is not shared by other ensembles. In 
other ensembles one can argue generally that such non-monotonic behavior cannot occur. 
For example, the conjugate distribution P{F) = J Q{r)e^''^df is monotonic. This follows 
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from noticing that the 3x3 matrix 

dlnP 



dFidFj 



=< in- < n >)(r,- < >) > (12) 



is positive definite. Such arguments do not apply to Q{f) since there is no analogous formula 
to Eq. (4) in the conjugate distribution. Indeed, if there were, a double humped form could 
not appear in Q(r), for one could express the second derivative of Q{r) as the variance of 
the force. Can P{z) show non-monotonic behavior? The answer is no, for it has been shown 
in [11] that 

— -1- = Qi^) 

z dz 

Since Q{r) is a probability density and therefore non-negative, it follows that ^ < for 
z > 0, thus ruling out multiple peaks in P{z). This argument which rules out multiple peaks 
in P{z) is a new observation. Note that a P{z) measurement difi:ers from Q{r) only in the 
transverse stiffness A of the trap. By tuning A we can permit fluctuations in the transverse 
direction and therefore destroy the phase transition present in the stiff A limit. One would 
expect to see a critical stiffness A = Ac below which the phase transition is destroyed. 
Alternately, one could tune the mean force and expect to see the phase transition vanishing 
below a critical mean force F — for a flxed value of P in the intermediate range of /3. We 
emphasize that the non-monotonic features of Qlf) in the semiflexible range (3 ~ 3.8 are 
predictions of the WLC model which can be tested against experiments. A single molecule 
with its ends confined in optical traps is expected to show this fiexible to rigid transition. 
The effect can be dramatic however, if a large number of molecules co-operatively show 
such a transition. One could attach the ends of a collection of semifiexible polymers to 
supramolecular layers [22] and detect the fiexible to rigid transition signalled by a change in 
the inter-layer separation via a suitable probe. It may be possible to exploit this dramatic 
transition from fiexible to rigid behavior in technological applications. 

If one considers the replica thermodynamic limit of a semifiexible polymer one sees that 
force-extension curves continue to remain distinct in the Hclmholtz and the Gibbs ensem- 
bles. So while interpreting a force-extension curve obtained for a collection of semifiexible 
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polymers suspended between two arrays of traps, one needs to know if the curve is obtained 
using the soft trap setup or a stiff trap setup. However, for /3 —>■ oo, which is the usual 
thermodynamic limit, the two ensembles give rise to the same force-extension curve. This 
observation is consistent with the fact that inequivalence of ensembles can survive at the 
thermodynamic limit for systems with long range interactions [7,8,23]. In the context of 
semiflexible polymers, the persistence length Lp plays the role of the range of interactions. 
There has been some work [24] on the thermodynamics of particle systems in the presence 
of external macroscopic fields in classical and quantal contexts. In these papers the authors 
have dealt with the macroscopic limit of the definition of pressure which is analogous to 
the thermodynamic limit of the definition of force in our work. In particular the authors of 
Ref. [24] discuss the connection between the values of the pressure defined by two different 
thermodynamic limit procedures: in the first, the system is confined successively in a se- 
quence of boxes which grows to fill up the whole space. In the second, the system is in an 
external potential similar to the present context. In the case of a semiflexible polymer the 
force corresponding to a given extension is the same in the Gibbs and Helmholtz ensembles 
only in the thermodynamic limit of /3 — > oo. This is analogous to the second procedure of 
taking the thermodynamic limit given in [24] , where they recover the thermodynamic notion 
of pressure from an underlying microscopic definition when they let the scale factor go to 
infinity in the macroscopic limit. The results of this paper are therefore, consistent with the 
general treatment in [24] . 

The fluctuation effects mentioned here also have some biological signiflcance. In partic- 
ular, the process of gene regulation involves interaction between parts of a DNA molecule 
which are about less than one persistence length apart 34nm) [21]. Over such short 
segments of the DNA fluctuation effects would be signiflcant. The replica thermodynamic 
limit would also play a role in the concrete biological context of a network of actin filaments 
forming the cytoskeletal structure. 
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